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ABSTRACT 



, Context. The hot X-ray emitting gas in clusters of galaxies is a very large repository of metals produced by supernovae. During the evolution 
• of clusters, billions of supernovae eject their material into this Intra-Cluster Medium (ICM). 
Aims. We aim to accurately measure the abundances in the ICM of many clusters and compare these data with metal yields produced by 
supernovae. With accurate abundances determined using this cluster sample we will be able to constrain supernova explosion mechanisms. 
t-H ' Methods. Using the data archive of the XMM-Newton X-ray observatory, we compile a sample of 22 clusters. We fit spectra extracted from 
, the core regions and determine the abundances of silicon, sulfur, argon, calcium, iron, and nickel. The abundances from the spectral fits are 

subsequently fitted to supernova yields determined from several supernova type la and core-collapse supernova models. 
' Results. We find that the argon and calcium abundances cannot be fitted with currently favoured supernova type la models. We obtain a major 
^ C* | improvement of the fit, when we use an empirically modified delayed-detonation model that is calibrated on the Tycho supernova remnant. The 
O |' two modified parameters are the density where the sound wave in the supernova turns into a shock and the ratio of the specific internal energies 
' of ions and electrons at the shock. Our fits also suggest that the core-collapse supernovae that contributed to the enrichment of the ICM had 
,~ , progenitors which were already enriched. 

t^J Conclusions. The Ar/Ca ratio in clusters is a good touchstone for determining the quality of type la models. The core-collapse contribution, 
r-j which is about 50% and not strongly dependent on the IMF or progenitor metallicity, does not have a significant impact on the Ar/Ca ratio. 
• • ■ The number ratio between supernova type la and core-collapse supernovae suggests that binary systems in the appropriate mass range are very 
efficient (~ 5-16%) in eventually forming supernova type la explosions. 

'x 

5^ ■ Key words. X-rays: galaxies: clusters - Galaxies: clusters: general - intergalactic medium - Galaxies: abundances - Supernovae: general - 
Nuclear reactions, nucleosynthesis, abundances 



1. Introduction 

Clusters of galaxies are the largest gravitationally bound ob- 
jects in the universe. About 80% of the baryonic matter in the 
clusters is in the form of hot X-ray emitting gas that has been 
continuously enriched with metals since the first massive stars 
exploded as supernovae. The abundances of elements in this 
hot Intra-Cluster Medium (ICM) therefore correspond to the 
time-integrated yield of the supernova products that reached 
the ICM. About 20-30% of the supernov a products is bein g 
locked up in stars in the member galaxies dLoewensteinl20 04 ) . 
Because of the huge mass of the accumulated metals in the 
ICM, clusters of galaxies provide a unique way to test nucle- 
osynthesis models of supernovae on a universal scale. 

Since the launch of the ASCA satellite, it has been possi- 
ble to do abundance studies using multiple elements. Several 



groups (e.g. |F inoguenov et al. 2000; Fukazawa et al. 2000; 
iFinoguenov et al.ll2.Q0lk IBaumgartner et al.ll2005h used ASCA 
observations of a sample of clusters to study the enrichment 
of the ICM. They were able to measure the abundances of 
iron, silicon, and sulfur. Also neon, argon, and calcium were 
sometimes detected, but with relatively low accuracy. The spa- 
tial distributions of iron and silicon indicated that the core of 
the clusters is dominated by SN la products (Fe), while the 
outer parts of the clusters appear to be dominated by core- 
collapse supernova products (O). Using ASCA observations 
some authors already tried to constrain the specific fl avour of 
supe rnova type la mode l s (e.g. Dupke & Arnaudl2001 ). Others, 



like IBaumgar tner et alj d2005l) used ASCA data to find that 



Population-Ill stars should play an important role i n the enrich 



ment of the ICM. However, this result is debated (Ide Plaa et al 
2006h . 
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With the XMM-Newton observatory (IJansen et al.l 120011) . 
which has both a better spectral resolution and a much larger 
effective area compared to ASCA, it is in principle possible 
to extend the number of detectable elements to nine. The first 
abundances determined from a sample of clusters observe d 
with XMM-Newton were published by iTamura et aD J2004h . 
The general picture from the ASCA samples was confirmed, 
except for the fact that the silicon and sulfur abundance show 
a centrally-peaked spatial profile like iron. The oxygen abun- 
dance appears to be more uniformly distributed in the clusters. 

Recently, IWerner et ail J2006I) and Ide Plaa et all d2006l) 
analysed deep XMM-Newton observations of the clusters 2A 
0335+096 and Sersic 159-03, respectively. They were able to 
accurately measure the global abundances of about nine ele- 
ments in the cluster and fit them using nucleosynthesis models 
for supernovae type la and core-collapse supernovae. The fits 
show that ~30% of all the supernovae in the cluster are type la 
and about 70% are core-collapse supernovae. In their data, they 
found a clear hint that the calcium abundance in these clusters 
is higher than expected. 



In this paper, we extend the approach of Ide Plaa et al 



(2006) and IWerner et al.l (120061) to a sample of 22 clusters ob- 
served with XMM-Newton. We aim to accurately measure the 
chemical abundances of all robustly detected elements and fit 
model yields of type la and core-collapse supernovae to the re- 
sults. Naturally, we discuss the anomalous calcium abundance. 

In our analysis we use Ho = 70 km s Mpc~', Q m = 0.3, 
and = 0.7. The elemental abundances presented in this 
paper are give n relative to the proto-solar abundances from 
lLoddersl(l2003l) . 



2. The sample and methodology 

We use XMM-Newton data from a sample of 22 clus- 
ters of galaxies in the redshift interval z=0-0.2. The clus- 
ters are selected primarily from the HIFLUGCS sample 
(Rei prich & Bohrin ger 2002), because this sample is well stud- 
ied and contains the brightest clusters in X-rays. We only select 
the observations with the best data quality. 



2.1. Sample selection 

Since its launch, XMM-Newton has been used to obtain more 
than 500 cluster and group observations. However, not all of 
these observations are suitable to use in a sensitive abundance 
study. We choose the following selection criteria to ensure that 
we select a clean and representative sample. 

- The redshift of the cluster is between z=0-0.2. We select 
only local clusters and assume that all of these clusters have 
a comparable enrichment history. 

- The cluster core fits inside the field-of-view of XMM- 
Newton. We need a region on the detector that is not heavily 
'polluted' with cluster emission to estimate the local back- 
ground. This excludes large extended nearby clusters like 
Virgo and Coma. 
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Fig. 1. Systematic effects in the temperature profile of the clus- 
ter of galaxies Sersic 159-03. We made the three different pro- 
files by subtracting three background spectra wit h different nor- 



malisation. The blank field backgrounds by iRead & Pon man 



(2003) were scaled up and down by a factor of 10% and were 
subtracted from the EPIC spectra. The impact to the tempera- 
ture profile is clear: the temperature is well determined in the 
bright core, but the background subtraction plays an important 
role in the outskirts. 

- The clusters have reported temperatures between ~2 and 
10 keV. We exclude groups of galaxies and extremely hot 
clusters. 

- The clusters are part o f the HIFLUGCS sample 
dReiprich & Bohringerll2002l) . 

- The observation must not suffer from a highly elevated level 
of soft-protons after flare removal. 

We found 22 clusters that meet these requirements. They 
are listed in Table Q] Together, they have a total exposure time 
of 690 ks. Roughly 40% of the datasets suffer from a high level 
of so-called residual soft-protons. Soft protons have energies 
comparable to X-ray photons. They can induce events in the 
detector that cannot be separated from X-ray induced events. 
When the soft-proton flux is high, they create a substantial ad- 
ditional background. The soft-proton flare filtering is in general 
not enough to correct for this. An elevated quiescent level is in 
some cases only detectable as a hard tail in the spectrum. We 
check the spectrum in an 8-1 1' annulus centred on the core of 
the cluster. If this spectrum shows an obvious hard tail (more 
than two times the model count rate at 10 keV), then we ex- 
clude the cluster from the sample. 

The properties of the clusters we selected are diverse. In 
Table Q] we list a few basic properties of the clusters in our 
sample. The redshifts lie in the range between z=0.0214 and 
z=0.184. Sixteen clusters contain a cooling core (cc). 

2.2. Methodology 

Data reductio n is done following a procedure that is extensively 
described in de Plaa et al. ( 20061) . We first reprocess the data 
with SAS version 6.5.0. Then, we filter out soft-proton flares 



J. de Plaa et al.: Constraining supernova models using the hot gas in clusters of galaxies 



3 



Table 1. Summary of the cluster properties of this sample. Data are taken from Table 3 and 4 in Reiprich & Bohringer (2002) 
apart from the classification. We list the following properties: (1) Heliocentric cluster redshift. (2) Column density of Galactic 
neutral hydrogen gas in units of 10 20 cm 2 . (3) ROSAT flux in the energy range 0.1-2.4 keV in units of 10~ n erg s~' cm~ 2 . (4) 
Luminosity in the energy range 0.1-2.4 keV in units of /i 5( 2 1 4 4 erg s" 1 . (5) X-ray temperat ure in keV(note: in this table 90% 
errors are used). (6) Cluster radius (R500) in Mpc as listed in lReiprich & Bohringerl ( l2002l) . (7) Extraction radius used in this 



analysis in arcmin (0.2/?5oo). (8) Effective XMM-Newton exposure time in ks. (9) Classification (cooling core [cc] or non-cooling 
core [non-cc]). 



Cluster 


z 


N„ 


fx 


L x 


kr 




^extr 


Exposure 


Class 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


2A 0335+096 


0.0349 


18.64 


9.16 


4.79 


3.01 +0.07 


1.15 ±0.02 


3.74 


114 


TTT. 

cc"-* 


A 85 


0.0556 


3.58 


7.43 


9.79 


6.9 + 0.4 


1.68 + 0.06 


3.42 


12 


cc"-* 


A 133 


0.0569 


1.60 


2.12 


2.94 


3.8_ 09 


1 74+0-30 
-0.16 


2.46 


20 


cc"- 6 


A 1651 


0.0860 


1.71 


2.54 


8.00 


6.1 ±0.4 


1.73 ±0.08 


2.26 


8 


cc"-* 


A 1689 


0.1840 


1.80 


1.45 


20.61 


9.2 + 0.3 


2.20 ± 0.06 


1.31 


36 


cc"-* 


A 1775 


0.0757 


1.00 


1.29 


3.18 


3 f.9+0.20 
•'•"^'-o.ii 


1.36 ±0.06 


4.04 


23 


non-cc a '' 


A 1795 


0.0616 


1.20 


6.27 


10.12 


7.8 + 1.0 


1.89 ±0.13 


3.46 


26 


CC oJ ' 


A 2029 


0.0767 


3.07 


6.94 


17.31 


9.1 ± 1.0 


2.01 ±0.12 


2.95 


11 


cc"-* 


A 2052 


0.0348 


2.90 


4.71 


2.45 


3.03 ± 0.04 


1.10 ±0.02 


3.59 


29 


cc"-* 


A 2199 


0.0302 


0.84 


10.64 


4.17 


4.10 + 0.08 


1.43 ±0.04 


5.38 


23 


cc"- 6 


A 2204 


0.1523 


5.94 


2.75 


26.94 


7.2 + 0.3 


1.82 ± 0.05 


1.32 


19 


cc"-* 


A 2589 


0.0416 


4.39 


2.59 


1.92 


3 T 22 


1 9Q+0.37 
1-i -0.22 


3.52 


22 


non-cc' 


A3112 


0.0750 


2.53 


3.10 


7.46 


5 3 +0J 


1.53 ±0.14 


2.29 


22 


cc a - h 


A 3530 


0.0544 


6.00 


0.99 


1.25 


3.9 + 0.3 


1.47 ±0.14 


3.06 


11 


non-cc 


A 3558 


0.0480 


3.63 


6.72 


6.62 


5.5+0.4 


1.55 ±0.07 


3.70 


43 


cc"-* 


A 3560 


0.0495 


3.92 


1.52 


1.60 


3.2 + 0.5 


1.14 ± 0.12 


2.60 


25 


non-cc rf 


A 3581 


0.0214 


4.26 


3.34 


0.66 


1.83 + 0.04 


0.87 ± 0.03 


4.63 


36 


cc"-* 


A 3827 


0.0980 


2.84 


1.96 


7.96 


7.1 ± 1.1 


9 9^+0.60 
-0.37 


2.57 


20 


non-cc 


A 3888 


0.1510 


1.20 


1.10 


10.51 


8.8 + 1.3 


2.5 ±0.3 


1.82 


23 


non-cc"'' 


A 4059 


0.0460 


1.10 


3.17 


2.87 


4.4 + 0.3 


1.40 ±0.06 


3.44 


14 


CC"'* 


MKW3S 


0.0450 


3.15 


3.30 


2.87 


3.7 + 0.2 


1.29 ±0.05 


3.24 


35 


cc"-* 


S 159-03 


0.0580 


1.85 


2.49 


3.60 


3 +L2 


1 99+O.23 


2.45 


113 


cc"-* 



" Peres et al. (1998), b White et al. ( 1997), c Buote & Lewis (2004), d Bardelli et al. (2002). 



which exceed the 2 sigma confidence level. The background 
subtraction is performed by subtracting a spectrum extracted 
from a closed filter observation that we scale to the instrumental 
noise level of the particular cluster observation. Cosmic back- 
ground components are included in the spectral fitting phase. 

We extract the spectra from a circular region around the 
core of the cluster. In order to sample comparable regions in 
all clusters, we choose a physical radius of O.2.R500. The values 
of Rsqq are taken from Reiprich & Bohringei ( 20021) . When we 
use the radii of 0.2/?5oo, we sample the dense core region of 
the cluster. The radii in arcmin for every cluster are listed in 
Table [U 

The spectral components of the background are fitted to a 
spectrum extracted from an 8-11' annulus. This region near the 
edge of the detector generally contains little cluster emission, 
because we select clusters with a relatively small angular size 
(see Sect. 12. U . This allows us to estimate the local background 
without a large bias due to cluster pollution. Small (10%) un- 
certainties in the background do not affect our analysis. Fig.Q] 
shows that the temperature can be robustly measured in the core 
of Sersic 159-03, even if the background estimate would be off 
by 10%. Therefore, we concentrate our analysis on the bright 
cluster cores. 



We fit heuristicly the background spectra with four com- 
ponents: two Collisional Ionisation Equilibrium (CIE) com- 



pone nts with temperatures 0.07 and 0.25 keV dde Plaa et al 



2006), a power-law component with r=1.41, and a normalisa- 
tion fixed to a flux valu e of 2.24 x 10" 11 erg cm -2 s" 1 deg~ 2 in 



the 2-10 keV range dDe Luca & Molendill2004l) . Another CIE 
component is added to fit any remaining cluster emission. The 
results of the background fits are then used in the fits to the 
spectra of the core. 

For the sp ectral fitting of the cluster spectra we use the 
wdem model (lKaastra et alJ l2004t Ide Plaa et alJ 120051) which 
proved to be most successful in fitting cluster cores (e.g . 
Kaastra et all |2004 IWerner et al] 120061 Ide Plaa et all l2006h . 
This model is a differential emission measure (DEM) model 
where the differential emission measure is distributed as a 
power law (dY/dT oc T 1 '") with a high (T mlix ) and low tem- 
perature cut-off (T„ in ). We fix T min to 0.1 times r max like in 



de Pl aa et alJ (j2006). We quote the emission-w eighted temper- 
ature kr mean of the distribution dde Plaa et alj|2006l) . 



3. Results 

In this section, we apply the wdem model to the EPIC spectra 
of the clusters in the sample. We fix the redshift to the value 
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Table 2. Basic properties of the sample of clusters obtained from a fit of MOS and pn data. Nu is given in units of 10 20 cm -2 , 
Y — J n e n H dV is in units of 10 66 cm" 3 , and kT max and &T mean are in keV. The a parameter is a measure for the slope of the 
emission-measure distribution. 



Cluster 


N H 


Y 


kr max 




a 


VT 

>^ L mean 


X 1 1 dof 


2A 0335 


25.71 ± 0.09 


18.76 ± 0.09 


3.486 ± 0.016 


0.360 


± 0.007 


7 7^7 4- 
Z. / J 1 ± 


u.uu 


2927 / 1300 


A 85 


3.25 ± 0.09 


30.1 +0.2 


6.80 ±0.18 


0.50 


±0.05 


^11-1- 


n 1 7 
u. i / 


976 / 739 


A 133 


1.72 ± 0.10 


9.25 + 0.10 


4.29 ± 0.08 


0.36 


±0.03 


^ An -i- 


U.Uo 


1339/775 


A 1651 


2.10 ± 0.16 


24.0 ± 0.4 


7.2 ± 0.9 


0.3 


±0.2 


^ Q -t- 
J.y + 


i n 


857 / 756 


A 1689 


1.96 ± 0.09 


88.4 ± 0.8 


13.0 ± 1.1 


0.52 


±0.19 


Q 7 -i- 
y. 1 + 


1 n 


1106/749 


A 1775 


0.48 + 0.12 


9.82 ± 0.12 


3.58 ± 0.19 


0.03 


±0.05 


j.j + 


n 7 
u.z 


1124/758 


A 1795 


1.08 ± 0.04 


40.92 ± 0.16 


7.05 ± 0.13 


0.55 


±0.04 


J.ZZ ± 


n 1 7 

U. 1Z 


1703 / 877 


A 2029 


3.23 ± 0.07 


72.6 + 0.5 


9.7 ± 0.4 


0.51 


±0.08 


7 3 _t_ 
/ . J ± 


U.J 


1279/759 


A 2052 


3.22 ± 0.06 


5.24 ± 0.04 


3.73 ± 0.04 


0.411 


±0.014 


7 SQ -i- 
Z.oy ± 


U.UJ 


1336/814 


A 2199 


1.16 ± 0.04 


11.06 ± 0.05 


4.93 ± 0.08 


0.32 


±0.03 


Q7 -i- 
J.y I + 


U.Uo 


3328 / 1922 


A 2204 


7.30 ± 0.13 


103.8 ± 0.9 


10.1 ±0.5 


1.5 


±0.4 


6.5 ± 


0.4 


982 / 778 


A 2589 


3.54 + 0.10 


5.067 ± 0.049 


3.5 ± 0.3 


0.00 


±0.07 


3.5 ± 


0.3 


1087/775 


A3112 


1.12 + 0.07 


24.57 ± 0.17 


5.76 ±0.11 


0.44 


±0.04 


4.42 ± 


0.11 


1221 /806 


A 3530 


6.1 + 0.3 


2.55 ± 0.06 


3.6 ± 0.4 


0.03 


±0.13 


3.6 ± 


0.6 


890 / 709 


A 3558 


3.5 + 0.07 


10.19 ± 0.06 


8.1 ±0.3 


0.61 


±0.11 


5.9 ± 


0.3 


1190/768 


A 3560 


3.0 + 0.2 


1.83 + 0.04 


3.3 ± 0.3 


0.00 


±0.04 


3.3 ± 


0.3 


996 / 729 


A 3581 


4.36 + 0.11 


1.85 ±0.02 


2.14 ± 0.02 


0.267 


± 0.008 


1.765 ± 


0.018 


1499 / 777 


A 3827 


2.36 + 0.12 


28.2 ±0.3 


8.0 ± 1.0 


0.25 


±0.19 


6.7 ± 


1.0 


1159/791 


A 3888 


0.43 + 0.12 


39.0 ± 0.5 


9.8 ± 1.7 


0.0 


±0.2 


9.8 ± 


2.6 


1018/777 


A 4059 


1.44 + 0.06 


7.85 ± 0.06 


4.33 ±0.19 


0.19 


±0.06 


3.7 ± 


0.2 


1820/ 1426 


MKW3s 


2.99 + 0.06 


8.35 ± 0.05 


4.35 ± 0.08 


0.30 


±0.03 


3.53 ± 


0.08 


1219/822 


S 159-03 


1.00 ± 0.03 


13.7 ± 0.06 


3.08 ± 0.02 


0.238 


±0.010 


2.59 ± 


0.02 


3072 / 1495 



given in iReiprich & Bohringerl ([2002) and leave the Galactic 
neutral hydrogen column density (/Vh) free in the fit. 



3.1. Basic properties of the sample 



Table |2] shows the fit results for Nh and the temperature struc- 
ture. The normalisation (Y), the maximum temperature of the 
DEM distribution (krmax), and the slope parameter a are di- 
rectly obtained from the wdem fit. The emission- weighted tem- 
perature (krmean), derived from the fit parameters, is a good 
indicator of the temperature of the cluster core. In our sam- 
ple these temperatures cover a range between roughly 1 .7 keV 
(A3581) and 9.8 keV (A3888). The sample is slightly biased to 
low-temperature clusters. 



3.2. Abundance determination 

From each fit to a cluster and from each instrument (MOS and 
pn), we obtain the elemental abundances of oxygen, neon, mag- 
nesium, silicon, sulfur, argon, calcium, iron, and nickel. All 
these abundances, however, can be subject to various system- 
atic effects. We know, for ex ample, that the oxygen and neon 



abund ances are problematic dde Plaa et al. 20061: Werner et al 



2006). The spectrum of the Galactic warm-hot X-ray emitting 
gas (e.g. the local hot bubble) also contains O vm lines that can- 
not be separately detected with the spectral resolution of EPIC. 
The brightest neon lines are blended with iron lines from the 
Fe-L complex near 1 keV, which makes an accurate determina- 
tion of the abundance difficult. Therefore, we do not use these 
two elements in the rest of our discussion. In the following sec- 
tions we discuss two other possible sources of systematic ef- 
fects and work towards a robust set of abundances. 



A fit with the wdem model does not always lead to an ac- 
ceptable x 2 - value (see Table |2]). This is largely due to system- 
atic errors between the MOS and pn detectors. We describe 
these systematic differences extensively in Sect. 13.21 A sec- 
ond reason for the high x 1 can be the complicated temperature 
structure that is often observed in cluster cores. Because the 
wdem model is just an empirical DEM model, the real temper- 
ature distribution in the core of the cluster may be somewhat 
different. Finally, because some weak lines are not yet in the 
atomic database, small positive residuals can arise in line-rich 



regio ns like, for example, the Fe-L complex (lBrickhouse et al 
20001) . 



3.2.1. MEKALvs. APEC 

A possible source of systematic effects is the plasma model 
that we use. T here is an alternativ e for the MEKAL-based code 
called APEC dSmifh et al.ll200lb . In Table [3] we show a com- 
parison of two spectral fits to the spectrum of Sersic 159-03 
using a two-temperature model. One fit is performed using the 
MEKAL based CIE model, and the other with APEC. 

In general, the differences between the MEKAL and APEC 
fits are minor. Both plasma models are equally well capa- 
ble of fitting the spectrum, which is clear from the x 2 values 
(3882/3091 [MEKAL] and 3874/3102 [APEC]). We note that 
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Table 3. Comparison of the fits using MEKAL and APEC. 
We fit two CIE models to the spectrum of Sersic 159-03. 
Normalisation (Y = f n e n//dV) is in units of 10 66 cirr 3 . 



Parameter 


MriJvAL 




y, 
1 1 


J.JXUJ 


5.13+0.11 


Y 2 


1.4±0.3 


2.20±0.07 


kr, 


1.99±0.06 


1.96±0.05 


kr 2 


4.4+0.4 


3.47±0.08 


Si 


0.315+0.013 


0.331±0.013 


s 


0.292+0.019 


0.310±0.019 


Ar 


0.20±0.05 


0.20±0.05 


Ca 


0.45±0.08 


0.51±0.05 


Fc 


0.504±0.007 


0.514±0.008 


Ni 


0.63+0.06 


0.71±0.05 


X 2 /dof 


3882/3091 


3874/3102 




1 .5 2 
Energy (keV) 



Fig. 2. Residuals showing the calibration difference between 
EPIC MOS and pn. We show the residuals of the cluster Sersic 
159-03 as a typical example. The data points in this plot show 
the difference between the pn data and the best-fit model of 
MOS: (observed[pn]-model[MOS])/model[MOS]. We indicate 
the (zero-redshift) line energies for the strongest affected lines 
in this band: Mg xn La ( 1 .47 keV), Si xn HE 4 ( 1 .87 keV) and 
SixivLa(2.01 keV). 



there is no standard routine to fit a wdem-like emission-measure 
distribution with APEC. Therefore, we test two-temperature 
models here, which is often also a good approximation for 
the temperature structure in a cluster core. The only differ- 
ences between the best-fit parameters of the two models are 
in the higher temperature component. MEKAL and APEC find 
a slightly different mix between the high and low temperature 
component. Despite the small differences in temperature struc- 
ture between the two codes, the derived abundances are con- 
sistent within errors. This concl usion is in line with a si milar 
comparison of the two codes bv lSanders & Fabian d2006l) . 



3.2.2. Cross-calibration issues 

EPIC cross-calibration efforts (Kirsch 2006) have shown that 
there are systematic differences in effective-area between MOS 
and pn that are of the order of 5-10% in certain bands. 
Systematic errors of this magnitude can have a large impact on 
abundance measurements. We fit the MOS and pn spectra sep- 
arately to investigate the impact on our abundance estimates. 

The main differences in calibration between MOS and pn 
can be found in the 0.3-2 keV band. In Fig. [2] we show an 
example of the differences we observe between the two instru- 
ments. The pn instrument shows a positive excess with respect 
to MOS in the 0.3 to 1.2 keV band. Note that the models fit 
well to the spectra from both instruments if the spectra are fit- 
ted separately. Between 1.2 to 2.2 keV, the pn gives a lower 
flux than MOS. These differences can have a significant effect 
on the abundances that are measured in this band, like for mag- 
nesium and silicon. 

We show the effect in Fig. fusing the spectrum of Sersic 
159-03. The plot shows the line contributions of the elements 
that contribute the most to the 1.0^+. 5 keV band. Between 
roughly 1 and 2.5 keV, the lines in pn appear to contain less 
flux than their equivalents in MOS, which is in line with the 
differences we found between the two instruments. The effect 
is most notable at 1 .4 keV, where the magnesium abundance is 
used by the fit to make up for the difference in flux. Because 
the flux of the iron feature at that energy is firmly coupled to 
the iron-K line complex, the magnesium line flux is the only 
one that can fill the gap in flux. The silicon line, however, is 
quite clean. But there is still a difference in flux at the position 
of this line. At higher energies between 2.5 and 4.5 keV, there 
is no significant effect anymore, which suggests that the sulfur 
and argon abundances are clean. Above 5 keV, there is a small 
difference in the slope that can influence the temperature and 
subsequently the calcium, iron, and nickel line fluxes. 

Because abundances are directly derived from these line 
fluxes, we should be able to see the differences in the mea- 
sured abundances. In Fig. |4] we show the abundance ratios for 
MOS and pn separately. We use the data of Sersic 159-03 as an 
example, because it is representative for the whole sample. The 
values for the sulfur, argon, and calcium abundances appear to 
be consistent within errors in the two instruments. However, 
the silicon and magnesium abundances are clearly not. From 
spectral fits to the pn spectra we obtain systematically lower 
abundances for silicon and magnesium relative to MOS. 

In order to check whether the systematic differences in the 
abundances are largely due to effective area effects, we fit the 
Sersic 159-03 spectra again with corrected effective areas. We 
correct the MOS effective area with a simple broad-band spline 
model such that it nearly matches the pn effective area over the 
whole band. We do the same for pn. The filled symbols in Fig. [4] 
show the corrected abundances. The corrected MOS abundance 
is consistent with the original pn abundance and vice versa. 
Therefore, we can conclude that the effective area is the main 
contributor to the systematic differences between abundances 
determined with MOS and pn. 

We choose to use a conservative approach and estimate the 
systematic error from Fig. |4] There are three elements that suf- 
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Fig. 3. Line contributions with respect to the best-fit continuum in the EPIC MOS (left) and pn (right) spectra of Sersic 159-03. 
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Fig. 4. Comparison between the abundances measured with the 
MOS and the pn instruments in the cluster Sersic 159-03. We 
plot the abundance ratios with respect to iron. The open sym- 
bols correspond to the original data. The filled symbols show 
the results when the effective area is corrected to match the ef- 
fective area of the other instrument. 

fer from systematic effects: magnesium, silicon, and nickel. 
The systematic error in magnesium has such a large magni- 
tude that we can not obtain a significant value for it. For silicon 
and nickel we calculate the weighted average and add the sys- 
tematic error linearly to the statistical error. The error should 
be large enough to cover both the MOS and pn results. Using 
this method, we derive the following systematic errors (relative 
error with respect to the average abundance value): Si (+11%) 
andNi(+19%). 

3.2.3. Intrinsic scatter 

In Fig. [5] we show the abundance ratios of S/Fe and Ca/Fe for 
the individual clusters. At first sight, the abundance ratios ap- 
pear to be consistent with being flat, but they have a small scat- 



ter. In principle we expect to find a scatter, because the clusters 
in our sample are morphologically different and may have had 
different chemical evolution history. The intrinsic differences 
between the clusters need to be taken into account if we can 
detect the scatter with high significance. 

In order to quantify this intrinsic scatter in the abundances, 
we calculate the error-weighted average of the abundances (see 
Table|l]i with an error as described in Eq.[T] 



222 
cr, n , = cr m + erf, 



0) 



The value for the measured uncertainty (cr m ) is known from the 
spectral fits, but the combined uncertainty (<x tot ) and the intrin- 
sic scatter in the population of clusters (o"i nt ) are yet to be deter- 
mined. We do this by varying <x int until the^ ed of the weighted 
average is equal to 1 . The variance in the x 2 distribution for n 
free parameters is 2n by definition. We use this variance to find 
the lcr limits on our estimate for <x int . The values we derive for 
(Tint are listed in Table [4] 

We find that the intrinsic scatter (cn nt ) in the silicon and 
sulfur abundance ratios differs significantly from zero. This in- 
trinsic scatter in the data needs to be included in the error of 
the weighted mean. Therefore, we use a new weighted mean 
for silicon and sulfur with -4- as weighing factor (see Tabled. 

°tot 

Presumably due to lower statistics the cr; nt of argon, calcium, 
and nickel does not show a significant deviation (> 3<x) from 
zero, hence we may employ the original weighted means. 

3.2.4. Final abundance ratios 

Now, we have derived values for the most relevant systematic 
uncertainties that affect our abundances. Using the total statis- 
tical uncertainty (cr tot ) and the uncertainty in the effective area, 
we calculate the final abundance values with their errors. This 
final set of abundance ratios is shown in Table [4] Silicon and 
nickel are both dominated by the systematic uncertainty in the 
effective area. The sulfur abundance is dominated by the intrin- 
sic scatter. 
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Fig. 5. Plot of the S-Fe ratio (left) and the Ca-Fe ratio (right). In the Ca/Fe plot we do not plot clusters with 2cr upper limits larger 
than 2.0 for plotting purposes (A1689, A3530, A3560, A3827 and A3888). 



Table 4. Wei ghted averages of the abundance ratios with re- 
spect to solar dLoddersll2003l) in our sample. Systematic uncer- 
tainties due to the effective-area calibration and uncertainties 
due to intrinsic scatter are included in the errors. The <x int that 
we list here is the intrinsic scatter per data point. 



X/Fe 


Weighted mean 


CTinl 




(incl. cr tot ) 




Si/Fe 


0.66±0.13 


0.17+0.05 f 


S/Fe 


0.60±0.06 


0.18±0.06 


Ar/Fe 


0.40±0.03 


0.11+0.05 


Ca/Fe 


1.03±0.04 


0.12+0.08 


Ni/Fe 


1.41+0.31 


0.2+0.2 f 



f o"i nt taken from PN data only. 



4. Discussion 

We determined the elemental abundances in the core of 22 clus- 
ters of galaxies with XMM-Newton. Most of the abundance s 
are not consistent with proto-solar abundances (Lodde rsl2003l) . 
The intrinsic scatter in the cluster abundance ratios is between 
0-30%, which is quite small. Our sample consists of both re- 
laxed and non-relaxed clusters as well as cooling and non- 
cooling core clusters. The small intrinsic scatter shows that 
the effects of merging, cooling and temperature structure on 
the abundance ratios is limited to 30% in cluster cores. We do 
not resolve a clear trend of abundances with the presence of a 
cooling-core. 

It is a well established idea that most of the metals from 
oxygen up to the iron group are generated by supernovae. We 
construct a few models using elemental yields of supernova 
type la (SNIa) and core-collapse supe rnova (SNcc) . This anal- 
ysis is simila r to the one described in IWerner et al. (2006) and 
dePlaa et alJd2006l) . 

We try several SNIa yields which we o btain from two phys - 
ically different sets of supernova models dlwamoto et al . 1999), 
namely slow deflagration and delayed detonation models. The 



W7 an W70 models describe a slow deflagration of the stellar 
core, while the other models are calculated using a delayed- 
detonation (DD) scenario. WDD2 is the currently favoured 
SNIa explosion scenario. 

For the core -collapse super novae we use the yields from a 



recent model bv lNomoto et a l. (2006). Note that with SNcc we 



mean all types of core-collapse supernovae including types II, 
lb, and Ic. We integrate the yields from the model over the stel- 
lar population using an Initial -Mass Function ( IMF). We per- 
form the calculation following lTsuiimoto et al.l d 1995b : 



Mi 



iioMg M,(m)m-('^dm 
J w „ u m-( 1+J£ ) dm 



(2) 







where Mfyri) is the /th element mass produced in a star of main- 
sequence mass m. We use a standard model with Salpeter IMF 
(x=1.35) and solar-metallicity (Z=0.02). 

For every element i the total number of atoms jV, is a lin- 
ear combination of the number of atoms produced by a single 
supernova type la (Y, j a ) and type cc (Y !jCC ). 



Ni = aY iM + bY u 



(3) 



where a and b are multiplicative factors of type SNIa and core- 
collapse supernovae respectively. The total number of particles 
for an element can be easily converted into a number abun- 
dance. This reduces to a system of two variables (a and b) and 
six data points (Si, S, Ar, Ca, Fe and Ni). The fits are inde- 
pendent of the values for the solar abundances, because they 
are divided out in the procedure. In essence, we fit the absolute 
abundances in the cluster. In the following sections, we present 
the ratio of the relative numbers of SNIa with respect to the 
total number of supernovae (SNIa + SNcc) that have enriched 
the ICM. 

The supernova number ratios that we present, reflect the 
supernova ratio that is fitted to the abundances that we mea- 
sure in the ICM. This does not necessarily correspond to the 
true ratio of supernova explosions in the entire cluster over 
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its lifetime dMatteucci & Chiappini 2005 ). Because SNIa ex- 
plode some time after the initial star burst, there may be a 
difference between the fractions of type la and core-collapse 
produ cts locked up into stars . However, this delay time (<3 



Gyr, iMaoz & Gal- Yam! 120041) is probably short with respect 



to the formation time scale of the cluster. Therefore, it is likely 
that the bulk of the metals formed before O.lfHubbie- Because 
this enrichment timescale is an order of magnitude smaller 
than the Hubble tim e, the instantaneous recycling approxima- 
tion (|Tinsley 1980) is presumably a reasonable approximation 
in the case of clusters. Using this approximation implies that 
we ignore the stellar lifetimes and thus the delay with which 
some chemical elements are released from stars into the ICM 



dMatteucci & Chiapp ini 2005). 

We can make a very rough estimate of the systematic uncer- 
tainty that we introduce in our supernova ratio when we adopt 
the instanteneous recycling approximation. Galactic evolution 
models provide an indication about how the O/F e ratio behaves 



in an d around galaxies (for example, Calura & Matteuccl 
2006). However, galactic models are based on specific assump- 
tions and approximations. Therefore, they also contain sys- 
tematic uncertainties that are not well known. The plots in 
Calura & Matteuccil (12006b . for example, suggest that the frac- 
tion of oxygen that is locked up in stars is about a factor of two 
higher then for iron. If this model is a reasonable representation 
of galactic evolution in clusters of galaxies, then we overesti- 
mate our SNIa/(SNIa+SNcc) ratio in clusters with respect to 
the true supernova ratio with about 40% at maximum. 



4.1. Solar abundances 

The abundance ratios of silicon, sulfur, and argon that we de- 
rive from the sampl e are lower the n proto-solar abundance ra- 
tios determined by (ILoddersll2003l) . If we fit a constant to the 
cluster abundance ratios, we obtain ax 2 of 418 for 5 degrees 
of freedom. This means that the chemical enrichment in cluster 
cores differs significantly from that in the solar neighbourhood. 

In order to compare the supernova ratios in clusters with the 
solar ratio, we fit the supernova models to a constructed dataset 
with solar abundance ratios that have a nominal error of 5% on 
every data point. From this fit (with a ^ 2 /dof of 64/4), we find 
a supernova type la (WDD2) contribution in the solar abun- 
dance of 0.15 (+ 0.08), which is actually simi lar to the value 
for our galaxy found by iTsuiimoto et aL ( 1995b . This suggests 
that the abundances in the Sun are probably dominated by core- 
collapse supernovae that usually produce nearly flat abundance 
ratios. However, by adopting the instantaneous recycling ap- 
proximation we might underestimate the type la fraction for the 
solar neighbourhood with the sam e factor that we derived f rom 
the galactic evolution models bv lCalura & Matteucci d2006l) . 



4.2. Supernova type la models by Iwamoto et al. 

We now try to fit the current supernova models to the data of the 
sample. In Table [5] we show the fit results using the supernova 
type la models W7, W70, WDD2, WDD3, and CDD2. None 
of the models provides an acceptable fit. The model with the 



Table 5. Number ratios of supernovae type la over the to- 
tal number of super novae der ived using SN la m odels by 
llwamoto et all d 1999b and by dBadenes et al J 12006). The re- 
sults are from a fit which contained the elements from silicon 
to nickel. We also list the results of the comparison to solar 
abundances. 



Model 


SNIa/SNIa+SNcc 




Constant 




418/5 


Solar 


0.15+0.08 


64/4 


W7 


0.22 ± 0.06 


152/4 


W70 


0.26 ± 0.07 


104/4 


WDD2 


0.37 ± 0.09 


84/4 


WDD3 


0.22 ± 0.06 


105/4 


CDD2 


0.32 + 0.08 


86/4 


Tycho 


0.72 + 0.17 


26/4 
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Fig. 6. Abundance ratios versus atomic numbers for the sample. 
We fit the supernova yield models for SN la (WDD2) and SNcc 
(Z=0.02 and Salpeter IMF). The black line shows the total fit, 
while the dotted and dashed lines represent the SN la and SNcc 
models respectively. 

lowest x 1 is the delayed-detonation model WDD2 (j^ 2 /dof = 
84/4). 

The reason why the models fail can be found in Fig. [6] 
The calcium abundance is highly underestimated by the mod- 
els. Moreover, the high calcium abundance forces the fit to in- 
crease the SNcc contribution. The current models are clearly 
not able to produce the observed Ar/Ca and Ca/Fe abundances. 
This result re-affirm s the earlier measure ments in 2A 0335+ 096 
dWerner et al.ll2006h and Sersic 159-03 dde Plaa et al.ll2006l) . 



4.3. Supernova type la models based on Tycho 



Recently, Badenes et al. d2006l) compared type la models by 



Chieffi & Stranierol d 1989b and iBravo et alj d 1996b to XMM 



Newton EPIC observations of the Tycho supernova remnant. 
The Tycho supernova is tho ught to have been a type la super- 
nova. iBadenes et al.l d2006J) empirically modified the parame- 
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Fig. 7. Same as Fig |6j but now w e fit the SNIa yields found 
in the Tycho supernova remnant bv lBadenes et alj (I2006I) . The 



nickel yield of the Tycho SNIa model was kindly provided by 
Carles Badenes (priv. comm.). 



ters of their delayed-detonation model to fit the Tycho obser- 
vations. The parameter that mainly determines the outcome of 
their model is the density where the subsonic wave, which runs 
through the white dwarf during the explosion, turns into a su- 
personic shock. This transition from deflagration to detonation 
is put in by hand in every DD model. By modifying this pa- 
rameter and the ratio between the specific internal energies of 
ions and electrons (J3), they found a best-fit delayed-detonation 
model that fitted the Tycho observations. 

We take the yields from this best-fit model of Tycho and use 
them as a supernova type la model in our fit to the cluster abun- 
dances. Note that in Tycho not all the ejected material is visible 
in X-rays, because the reverse shock has not ionised all the ma- 
terial yet. Therefore, the Tycho results might not reflect the total 
SNIa yields yet. Despite this caveat, the Tycho model provides 
a majo r improvement in x 2 compared to the llwamoto et al 



( 1999) models (see Table©. In Fig. El we show that the Tycho 
model is more successful in fitting the calcium abundance in 
clusters. Moreover, the supernova ratios change dramatically, 
the SNIa/(SNIa+SNcc) ratio for this model is 0.72 ±0.17. 

This shows that the Ar/Ca and Ca/Fe abundance ratios 
mainly determine how well type la models fit. By varying the 
parameters of the delayed-detonation models, it is in principle 
possible to obtain a calcium abundance that fits the observa- 
tions, which apparently very effectively constrain type la mod- 
els. 

4.4. Core-collapse models 

In order to test whether our models are also reproducing the 
core-collapse contribution adequately, we need abundances 
of some typical core-collapse products. Therefore, we esti- 
mate the oxygen and neon abundance of the sample using the 
Reflection Grating Spectrometer (RGS) aboard XMM-Newton. 
For oxygen, we take the average of the clusters Sersic 159-03 
and 2A 0335+096, because these clusters have the highest ex- 
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Fig. 8. Fit using the SNIa yields by iBadenes et al. 
but now with additional oxygen and neon data points (stars) 
obtained from the RGS spectra of Sersic 159-03 and 2A 
0335+096. Here, the core-collapse model with Z=0.02 and 
Salpeter IMF is used. 



posure in our sample an d very good RGS data dde Plaa et al. 
20061: Wern er et al . 2006). The O/Fe measurements of the two 
clusters are not statistically consistent possibly due to sys tem- 
atic differences in the line widths (see lde Plaa et al. 20061 for 
an explanation of this effect). Therefore, we take the average 
value and assign an error which covers both results within lcr. 
The neon abundance is consistent in both Sersic 159-03 and 
2A 0335+096, hence we take the weighted average of the two 
neon abundances and use them in the rest of the fits. 

The fit including O/Fe and Ne/Fe from RGS is shown in 
Fig. [8] We use here a standard core-collapse model (Z=0.02 
and Salpeter IMF) and the type la model based on Tycho. The 
trend in O/Ne that is predicted by the core-collapse model, is 
consistent with the O/Ne ratio that we observe. However, the 
core-collapse contribution needs to increase with respect to the 
model used in Fig.Efro explain the absolute values for O/Fe and 
Ne/Fe (see Tabled- The typical values that we derive are of the 
order of 45-60%. Considering the uncertainties, this number 
is compatible with the current supernova type la ratio within 
z=0.03 (~42%) determined by the Lick Obs ervatory Supernova 
Search (LOSS) (Ivan den Bergh et alj|2005h . 

In principle, the increase of the core-collapse contribution 
with respect to the results before we included oxygen an neon 
results in a smaller predicted Ar/Ca ratio. However, the data 
suggest that the Ar/Ca ratio is larger. The plot shows that 
this particular core-collapse model still allows a relatively high 
Ar/Ca ratio, because the absolute contributions of silicon, sul- 
fur, argon, and calcium are relatively small in this model. 

4.4.1 . Effect of progenitor metallicity on core-collapse 
yields 

Up to now, we have used a core-collapse supernova model 
that assumes that the progenitor had a solar metallicity. 
Nomo to et al. (2006) also provide models where the metallic- 
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Fig. 9. Fit using the SNIa yields by iBadenes et aL ( 2006b and Fig. 10. Fit using the SNIa yields by Badenes et al.1 (2006) and 
SNcc yields for Z=0. 001 metallicitv progenitors ( Nomoto et al. SNcc yields for Z=0. 001 metallicitv progenitors ( Nomoto et all 
2006). 120061) that are integrated over a "top-heavy" IMF 



Table 6. Results for a su pernova fit to the a bundance ra- 
tios using SNcc models by iNomoto et alj (12006b that are in- 
tegrated over a Salpeter (jt=1.35) or "top-heavy" IMF (x=0). 
Z=0.02 is the solar abundance. For the type la model we used 
IBadenes et all J2006I) . 



Salpeter IMF 



"Top-heavy" IMF 





SNIa 


x'/dof 


SNIa 


*7dof 


SNIa+SNcc 


SNIa+SNcc 





0.55+0.05 


79/6 


0.68±0.04 


102/6 


0.001 


0.51+0.05 


54/6 


0.62±0.04 


69/6 


0.004 


0.49±0.05 


55/6 


0.62±0.04 


58/6 


0.02 


0.44±0.05 


40/6 


0.57±0.04 


34/6 



ity (Z) of the supernova progenitor is 0, 0.001, and 0.004. In 
Table [6] we show the results for the fits using progenitor metal - 
licities ranging from to 0.02 (solar). The Z=0.001, Z=0.004, 
and Z=0.02 models do show a relatively small variation in^ 2 . 
The data are still compatible with a wide range of metallicities 
(0.001-0.02). 

In Fig.|9]we show the fit result for the Z=0.001 model. The 
main difference between this Z=0.001 model and the Z=0.02 
model is the amount of oxygen and neon produced. The neon 
peak is clearly less pronounced compared to the model in 
Fig. [8] while the plateau from silicon to nickel in the core- 
collapse contribution is nearly unaffected. That also confirms 
that the influence of metallicity differences in SNcc models on 
the Ar/Ca ratio is limited. 

4.4.2. Effect of IMF on core-collapse yields 

We also fit the data for core-collapse models integrated over a 
"top-heavy" IMF (x=0), because presumably more high-mass 
stars form in low-metallicity environments. Table [6] shows that 
the fits are similar to the fits using Salpeter IMF. In Fig. \W\ 
we show the fit for this metallicity. The main differences with 



Salpeter IMF models are in the O/Ne ratio. The abundance peak 
in the model is shifted from neon to oxygen with respect to the 
Salpeter models, which is less consistent with the neon abun- 
dance from RGS. Again, the plateau from silicon to nickel in 
the core-collapse model is barely affected. 

4.5. The fraction of low mass stars that become Type 
la SNe 

This study indicates a much higher lifetime averaged ratio of 
type la to core-collapse supernovae than in the Galaxy. The rea- 
son for this large contribution of SNIa is likely that in late type 
galaxies, like our own galaxy, ongoing star formation ensures 
an ongoing core collapse contribution. For clusters of galaxies, 
the star formation continued at a very reduced level shortly af- 
ter the formation of the cluster. This has some interesting con- 
sequences. SNIa are likely the result of a thermonuclear run- 
away explosion of a C-O white dwarf in binary systems, caused 
by accretion from the secondary star. Since this involves the 
formation of a white dwarf from a star with mass M < 10 Mq, 
there is a considerable delay between the period of star forma- 
tion and the resulting SNIa explosion. In late type galaxies the 
subsequent waves of star formation make it difficult to disen- 
tangle the SNIa contribution from recent and old star formation 
periods. However, cluster of galaxies are an interesting labora- 
tory to study the fraction of all stars that will eventually become 
type la supernovae, since the star formation has shut down, and 
for the last few Gyr since formation, the buffer of potential type 
la progenitors has nearly emptied. 

As a result, the fraction of SNIa in clusters must be a good 
approximation to the fraction of low mass stars that can become 
type la ex plosions. For a po wer law initial mass function we 



type la explosions, ror a pov 
can write (lYoshii et al.l [T996): 



SNIa 



SNIa + SN CC fsma j*f<* m -d+x) dm + m -a+x)dm 



(4) 
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Table 7. The fraction, fsNia, of low mass stars in the range 
M/ ow to M cc that will eventually result in SNIa. It is based on 



our observed SNIa fraction of 



SNIa+SNcc 



= 0.44 + 0.10. 



M w = 0.9M o M low = 1.5 M© 



Salpeter IMF (x = 1.35) 
M cc = 8 M© 4.1 ± 1.7% 9 ± 4% 

M CC = 10M© 2.9 ±0.1% 6 ±3% 

Kroupa (2002) 
M cc = 8 M© 2.0 ± 0.8% 4.0 ±1.7% 

M cc = 10 M© 1 .3 ± 0.6% 2.7 ± 1.1 % 



with M) ow the lower limit to the mass of stars that can have con- 
tributed to the type la production, and M up the mass of the most 
massive stars. The parameter of interest here is fsNia, which is 
defined as the fraction of stars with M\ ow < M < M cc that will 
eventually explode as SNIa. 

It is clear that the absolute lower limit to M/ oll , is the 
mass of stars that can evolve to a C-O white dwarf within a 
Hubble time, about 0.9 Mq. However, in order to have suf- 
ficient mass available in the binary to push the white dwarf 
over the Chandrasekhar limit, it is usually assumed that the to- 
tal mass of a binary producing a SNI a explosion should exceed 
~ 3 M© (Mat teucci & RecchfcoOlh . This implies that the ini- 
tial mass of the primary star > 1.5 Mq. 

The time scale for producing a Type la supernova is 
most likely determined by the evolution of the secondary, i.e. 
roughly by the duration of the main sequence. If the secondary 
is 0.9 Mq, this corresponds roughly to a Hubble time, but there 
is considerable evidence that the mean delay time between star 
formation and type la explo s ion is shorter, namely of the or- 
der of 1 Gyr. IStrolger et all d2004 find a value of 2-4 Gyr. 
How ever, based on a lack of o bserved SNIa in clusters of galax- 
ies, 



Maoz & Gal-YarrJ |2004 find a 2< x upper-limit of a bout 3 

<l2006h 



Gyr for the delay time. Very recently Ma nnucci et al 



argued for two channels for SNIa explosions, one with a very 
short delay time of 10 8 yr, and the other with 2-4 Gyr. In other 
words, the majority of the type la explosions must occur in bi- 
naries where the secondary star has a lifetime shorter than 2-4 
Gyr. The mass of the secondary star that corresponds to such 
lifetimes is about 1.25-1.5 M0 (IGreggio & Renzinilll983l) . In 
the extreme case, where the delay time is < 2 Gyr, M/ OV) would 
be about 1.5 M0. 

Table [7] shows the values for fsNia using different assump- 
tions for M/ ow , M cc and the initial mass function. We use the 
best fit SNIa rate of 0.44±0. 10 (Tabled, but with slightly larger 
errors in order to allow for systematic uncertainties. The values 
that we find for M/ OM = 1.5Mq are high compared to previ- 
ous results based on other observational data, which suggests 
f SNIn = 1- 5% (lGreggiq&Renzini|[l9^lYoshiiet^ri996l; 



Matteucci & RecchfeOOlh . Note that these authors sometimes 
use slightly different definitions for M\ ow and M cc . Our low- 
est value is 1.3%, which assumes M cc = 10 M0 and Mi„ w = 
0.9 M (T) together with the broken power law IMF of iRroupa 
(120021) . Notice also that fsNia refers to all stars with masses 



between M\ ow and M cc . We can estimate the probability that 
a binary produces a type la supernova from the value of f$ ma ■ 
We assume that roughly 50% of all stars are formed in binaries, 
which introduces a factor of two. In addition, we need another 
factor of two because we want to count binaries instead of stars. 
Together, the probability that a binary produces a SNIa super- 
nova is thus about four times higher then fsNia- 

Our derived supernova ratios therefore suggest that binary 
systems in the appropriate mass range are very efficient in even- 
tually forming SNIa explosions (~ 5-16%, depending on the 
assumptions for the IMF and M cc ). We are aware that we ig- 
nore several complications in this simple calculation, such as 
an increased binarity fraction for massive stars that flattens the 
IMF for binary stars and binary mass ratios that may peak near 
1 . Also the instanteneous recycling approximation may intro- 
duce an additional uncertainty. However, a detailed calculation 
is beyond the scope of this paper. 

5. Conclusions 

We measure the abundances for silicon, sulfur, argon, calcium, 
iron, and nickel in a sample of clusters with XMM-Newton 
(EPIC), and we add a high-resolution oxygen and neon mea - 
surement from RGS dde Plaa et al]|200d [Werner et alj|2006l) . 
From these data we conclude that: 

- The Ar/Ca ratio in clusters is a good touchstone for de- 
termining the quality of type la models. The core-collapse 
contribution, which is about 50% and not strongly depen- 
dent on the IMF or progenitor metallicity, does not have a 
significant impact on the Ar/Ca ratio 

- Current supernova type la models dlwamoto et al .1 19991) do 
not agree with our data, because they fail to produce the 
Ar/Ca and Ca/Fe abundance ratios. 

- A major improvement of the supernova fits is obtained, 
when we use an empirically-modified supernova type la 
mode l, which is calibrat ed on the Tycho supernova rem- 
nant (Ba denes et al. I l2006h . This model largely solves the 
problems with the Ar/Ca and Ca/Fe abundance ratios by 
varying the density where the sound wave in the supernova 
turns into a shock and varying the ratio of the specific in- 
ternal energies of ions and electrons at the shock. 

- The number ratio between supernova type la and core- 
collapse supernovae suggests that binary systems in the ap- 
propriate mass range are very efficient (~ 5-16%) in even- 
tually forming supernova type la explosions. 

- We find that the progenitors of the core-collapse super- 
novae which contributed to the ICM abundances have prob- 
ably been enriched. Progenitor abundances range from 
Z -0.001 to Z -0.02. 

- The intrinsic spread in abundance ratios between clusters is 
smaller than 30%. That means that the chemical histories 
of the clusters do not depend a lot on cluster temperature, 
temperature structure or merging activity. 
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Appendix A: Abundance data 



In Table I A. 1 1 we list the abundances obtained from fits to the 
EPIC data. The MOS and pn spectra are fitted simultaneously. 
We included a systematic error in the uncertainties on the Si/Fe 
and Ni/Fe abundance ratios (see Sect. 13.2.21 for a discussion 
about systematic errors). 

List of Objects 

'2A 0335+096' on page 
'Sersic 159-03' on page 
'2A 0335+096'_on page 
'A 85' on page 
'A 133' on pag 
'A 1651' on page 
'A 1689' on page 
'A 1775' on page 
'A 1795' on page 
'A 2029' on page 
'A 2052' on page 
'A 2199' on page 
pdA 2204' on page 
'A 2589' on page 
'A 31 12' on page 
'A 3530' on page 
'A 3558' on page 
'A 3560' on page 
'A 3581' on page 
'A 3827' on page 
'A 3888' on page 
'A 4059' on page|3 
'MKW3S' on page]3j 
'S 159-03' on page^ 
'Tycho' on page|8j 
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Table A.l. Abundance ratios in the sample of clusters with respect to the solar abundances determined bv lLodd ers (20031). The 
listed ratio is calculated using F ^ F ^ 
discussion about systematic errors. 



We included a systematic error in the data points of Si/Fe and Ni/Fe. See Sect. l3.2.2l for a 



Cluster 


Si/Fe 


S/Fe 


Ar/Fe 


Ca/Fe 


Ni/Fe 


Fe 


2A 0335 


0.78 ± 0.09 


0.636 


± 0.019 


0.43 


± 0.04 


0.95 ± 0.06 


1.4 ± 


0.4 


0.741 


± 0.006 


A 85 


0.72 + 0.18 


0.61 


± 0.15 


0.4 


± 0.4 


1.4 ± 0.4 


1.0 ± 


0.7 


0.574 


±0.018 


A 133 


0.64 ±0.14 


0.40 


±0.09 


0.6 


± 0.2 


1.3 ±0.3 


1.6 ± 


0.6 


0.81 


±0.02 


A 1651 


0.0 ± 0.4 


0.3 


± 0.4 


0.0 


± 0.4 


0.0 ± 0.6 


1.4 ± 


1.3 


0.45 


±0.03 


A 1689 


0.3 ± 0.4 


1.1 


± 0.4 


0.8 


± 1.0 


0.6 ± 0.9 


0.6 ± 


1.0 


0.40 


±0.02 


A 1775 


0.57 + 0.18 


0.77 


± 0.14 


0.5 


±0.3 


1.7 ± 0.4 


1.6 ± 


0.7 


0.63 


±0.02 


A 1795 


0.75+0.14 


0.34 


±0.08 


0.2 


± 0.2 


1.1 ± 0.3 


1.3 ± 


0.5 


0.517 


± 0.009 


A 2029 


0.4 + 0.2 


0.32 


± 0.17 


0.00 


±0.10 


1.3 ±0.5 


1.6 ± 


0.7 


0.587 


± 0.017 


A 2052 


0.74 + 0.12 


0.76 


±0.05 


0.55 


± 0.11 


1.35 ± 0.16 


1.4 ± 


0.4 


0.682 


±0.011 


A 2199 


0.73+0.13 


0.60 


± 0.07 


0.60 


±0.16 


1.04 ± 0.19 


1.5 ± 


0.5 


0.532 


± 0.008 


A 2204 


0.75+0.18 


1.32 


± 0.19 


0.0 


± 0.2 


1.6 ±0.6 


1.6 ± 


0.7 


0.59 


±0.02 


A 2589 


0.55 ±0.14 


0.58 


±0.10 


0.4 


± 0.2 


1.2 ± 0.3 


1.8 ± 


0.6 


0.666 


±0.019 


A 3112 


0.70 ±0.14 


0.67 


±0.09 


0.5 


± 0.2 


1.3 ± 0.3 


1.9 ± 


0.6 


0.695 


±0.016 


A 3530 


1.1 ±0.6 


1.1 


±0.6 


0.3 


±0.9 


1.2 ± 1.5 


0.0 ± 


2.0 


0.28 


±0.04 


A 3558 


0.74 ±0.16 


0.41 


± 0.12 


0.3 


±0.3 


1.6 ± 0.4 


1.5 ± 


0.6 


0.478 


± 0.017 


A 3560 


0.9 ± 0.3 


0.4 


±0.3 


0.2 


±0.5 


1.0 ± 1.0 


1.1 ± 


1.3 


0.39 


±0.03 


A 3581 


0.80 ±0.10 


0.75 


± 0.04 


0.66 


±0.09 


1.32 ± 0.16 


1.3 ± 


0.4 


0.654 


±0.013 


A 3827 


0.2 ± 0.3 


0.0 


± 0.4 


0.0 


±0.7 


2.9 ± 1.2 


1.1 ± 


1.2 


0.38 


±0.03 


A 3888 


0.1 ±0.5 


1.8 


± 1.1 


0.0 


± 1.3 


1.0 ± 1.9 


0.0 ± 


1.8 


0.30 


±0.03 


A 4059 


0.59 ±0.13 


0.45 


± 0.07 


0.16 


± 0.17 


0.7 ± 0.2 


1.0 ± 


0.5 


0.687 


± 0.014 


MKW3s 


0.84 ±0.13 


0.63 


± 0.07 


0.54 


±0.18 


1.2 ± 0.2 


1.4 ± 


0.5 


0.551 


± 0.011 


S 159-03 


0.67 ±0.10 


0.55 


± 0.03 


0.41 


±0.07 


0.87 ± 0.10 


1.6 ± 


0.4 


0.533 


± 0.005 



